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Abstract 

Using Painleve analysis, the Hirota multi-linear method and a direct ansatz tech- 
nique, we study analytic solutions of the (l+l)-dimensional complex cubic and 
quintic Swift-Hohenberg equations. We consider both standard and generalized ver- 
sions of these equations. We have found that a number of exact solutions exist to 
each of these equations, provided that the coefficients are constrained by certain re- 
lations. The set of solutions include particular types of solitary wave solutions, hole 
(dark soliton) solutions and periodic solutions in terms of elliptic Jacobi functions 
and the Weierstrass p function. Although these solutions represent only a small 
subset of the large variety of possible solutions admitted by the complex cubic and 
quintic Swift-Hohenberg equations, those presented here are the first examples of 
exact analytic solutions found thus far. 
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1 Introduction 



Complicated pattern-forming dissipative systems can be described by the 
Swift-Hohenberg (S-H) equation [1,2]. Classic examples are the Rayleigh- 
Benard problem of convection in a horizontal fluid layer in the gravitational 
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field [3,4], Taylor-Couette flow [5], some chemical reactions [6] and large-scale 
flows and spiral core instabilities [7] . Examples in optics include synchronously- 
pumped optical parametric oscillators [8], three- level broad-area cascade lasers 
[9] and large aspect-ratio lasers [10,11,12,13]. In the context of large-aperture 
lasers with small detuning between the atomic and cavity frequencies, the 
complex cubic Swift-Hohenberg (CCSH) equation was derived asymptotically 
[14,15,16]. This equation is also believed to be relevant for oscillatory convec- 
tion in binary fluids. The appearance of spatial patterns is the most remarkable 
feature of the solutions. 

Despite their visual complexity, spatial patterns are actually formed from a 
certain number of ordered or chaotic combinations in space of some simpler 
localized structures. It is known that the S-H equation admits the existence of 
transverse localized structures and phase domains [14]. These coherent struc- 
tures can be considered as bright or dark solitons. Consequently, these simple 
localized structures and their stability are of great interest in the study of 
any pattern-forming system. If we keep these visual properties of patterns in 
mind, we can split the problem into several steps. Firstly, we should study 
the simplest stationary objects (localized structures). Secondly, we have to 
investigate their stability. As a last step, we should study their interaction 
and the combination structures. The distinction between stationary and mov- 
ing patterns is defined by the nature of the interaction forces between the 
localized structures. The above program cannot be carried out in one move. 
Moreover, in studying stationary localized structures, we have to start with 
the simplest cases, viz. (l+l)-dimensional structures. This sequential step-by- 
step approach allows us to avoid any possible confusion inherent in trying to 
explain a complicated structure in one step. 

Whichever part of the above program we consider, it is clear that it can mainly 
be done only using computer simulations. This has been done in the majority 
of papers published so far. We will analyse only (l+l)-dimensional bright and 
dark solitons, and leave aside, for the moment, the question of their stabil- 
ity and the interaction between them. Moreover, our task here is to obtain 
some analytic results. The latter point will further restrict the class of solu- 
tions which are of interest in our investigation. Nonetheless, this is a step of 
paramount importance because, until now, no analytic solutions have been 
known. 

In many respects, localized structures of the S-H equation are similar to those 
observed in systems described by the previously-studied complex Ginzburg- 
Landau equation(CGLE). We recall that some analytic solutions for the cubic 
and quintic CGLE are known [17,18,19,20]. For the cubic (1 + 1)-D equation, 
analytic solutions describe all possible bright and dark soliton solutions. In 
contrast, in the case of the quintic equation, the analytic solutions of the 
CGLE represent only a small subclass of its soliton solutions. Moreover, the 
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stable soliton solutions are located outside of this subclass [20]. Therefore, 
practically useful results can only be obtained numerically. 

We note that, apart from some exceptions, the CGLE generally has only iso- 
lated solutions [20,21] i.e. they are fixed for any particular set of equation 
parameters. This property is fundamental for the whole set of localized so- 
lutions of the CGLE. The existence of isolated solutions is one of the basic 
features of dissipative systems in general. The qualitative physical foundations 
of this property are given in [22] . This property is one of the reasons why the 
above program is possible at all. Like the CGLE, the S-H equation models 
dissipative systems, and we expect that it will have this property. Indeed, 
preliminary numerical simulations support this conjecture. 

The main difference between the S-H equation [23] and the CGLE lies in its 
more involved diffraction term. The latter is important in describing more 
detailed features of an actual physical problem. However, these complications 
prevent us from easily analyzing the solutions. In fact, it was not clear that 
such solutions could exist at all [8,11]. In this work, we study the quintic 
complex S-H equation in ID and report various new exact solutions. 

Before going into further details, we should distinguish between the real Swift- 
Hohenberg(RSH) equation and the CCSH equation. The former can be con- 
sidered as a particular case of the latter. Moreover, the RSH equation is a 
phenomenological model which cannot be rigorously derived from the origi- 
nal equations [3]. In contrast, the CCSH equation is derived asymptotically 
and is rigorous in the long-wavelength limit. In standard notations, the CCSH 
equation can be written as 

i(j t = rip - (1 + ic)|^| 2 ^ + iaAif) - d(A + A) 2 ip , (1) 

where r is the control parameter and a characterizes the diffraction properties 
of the active medium. In the limit of A — > 0, i.e., in the long- wavelength limit, 
the differential nonlinearities formally have higher orders and therefore can be 
dropped. In this case the wave-number-selecting term, (A + A) 2 i{j, is just a 
small correction to the diffraction term iaAip. Thus the CCSH equation can 
be treated as a perturbed cubic CGLE. For us the main interest here is (1 + 1) 
dimensional case when A has only x-derivatives. 

In physical problems, the quintic nonlinearity can be of equal or even higher 
importance to the cubic one [24] as it is responsible for stability of localized 
solutions. Sakaguchi and Brand made a numerical investigation of the complex 
quintic Swift-Hohenberg (CQSH) equation 

i> t = ai> + b\ifj\ 2 ij - c#| V - d(l + d xx fi) + ifd xx iP , (2) 
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where all coefficients a = oi + ia 2 , b = b\ + ib 2 , c = c\ + ic 2 , and the order 
parameter ifj are complex, but d and / are real [23]. For a, b and d real and 
c = and / = 0, the original Swift-Hohenberg equation, which was derived 
as an order parameter equation for the onset of Rayleigh-Benard convection 
in a simple fluid, is recovered. Then eq. (2) can be rewritten as 

iipt + (/ + 2id)i) xx + idif) xxxx + (b 2 - ibi)\tl)\ 2 ip + (-c 2 + zci)|V>| 4 V> 

— (— a 2 + i(a± — d))ip . (3) 

This equation can be generalized to 

%i>t + (3ip xx + Jipxxxx + Atl^l 2 ^ + ^IV'IV = iSip , (4) 



where all coefficients (3 — (3± + i(3 2 , 7 = 71 + ^72 , \i = A*i + ^2 and z/ = 27 + zz/ 2 
are complex and 5 is real. Eq.(4) is written in the form which we will refer 
to as the generalized CQSH (or GCQSH) equation in the rest of this study. 
Then the CCHS, real quintic Swift-Hohenberg (RQSH) and real cubic Swift- 
Hohenberg (RCSH) equations are particular cases of eq.(4) with some of the 
coefficients being equal to zero. 

The interpretation of the variables depends on the particular problem. In 
optics, t is the propagation distance or the cavity round-trip number(treated 
as a continuous variable), x is the transverse variable, (3\ is the 2nd order 
diffraction, is the 4th order diffraction, \x 2 is a nonlinear gain (or 2-photon 
absorption if negative) and 5 represents a difference between linear gain and 
loss. The angular spectral gain is represented by the coefficients /3 2 and 72. 

Sakaguchi and Brand observed some soliton-like structures in their numerical 
work [23]. They also showed that the CQSH equation admits a stable hole 
solution. Their results indicate that exact solutions may exist, but till now 
none have been found. In this paper, we study the CQSH and CCSH equations 
by using Painleve analysis and the Hirota multi-linear method. Our method 
is an extension of the modified Hirota method used earlier by Nozaki and 
Bekki [25] and the modification of the Berloff- Howard method [26]. We note 
that Nozaki and Bekki solved the CGLE by using the Hirota bilinear method, 
while Berloff and Howard solved some real-coefficient non-integrable equations 
using the Weiss- Tabor- Carnevale (WTC) method and the Hirota multi-linear 
method. In addition, we confirm that our solutions are valid by using a direct 
ansatz approach. The solutions we obtain here can be considered as a basis 
for further development. 

The rest of the paper is organized as follows. The methodology and analytical 
procedures are described in Sec. 2. Painleve analysis of the CQSH and CCSH 
equations is performed in Sec. 3. Exact solutions of the generalized CQSH and 



4 



CCSH equations and particular cases of them are obtained in Sec. 4. Finally, 
we summarize with our conclusions in Sec. 5. 



2 Methodology 

2.1 Painleve analysis and Hirota multi-linear method 

Any solution of an equation must be in accord with the singularity structure of 
that equation. Painleve analysis [27] is a tool for investigating that structure. 
This analysis can be applied both for ordinary differential equations(ODEs) 
and partial differential equations(PDEs). The power of the Painleve test lies 
in its easy algorithmic implementability. The main requirement is the repre- 
sentation of any possible solution in the form of a Laurent expansion in the 
neighborhood of a movable singularity: 



where r is the leading-order exponent, F(= 0) is a singularity manifold given 
by F(z), and Uj is a set of analytic functions of z. 

There are two necessary conditions for an ODE to pass the Painleve test: 

• the leading-order r must be an integer, and 

• it must be possible to solve the recursion relation for the coefficients Uj(z) 
consistently to any order. 

The general expansion of a non-integrable equation will fail the Painleve test 
at one of these two steps. Leading-order analysis can be done by balancing 
the highest-order derivative in x with the strongest nonlinearity. 

Weiss et al. developed the singular manifold method to introduce the Painleve 
property into the theory of PDEs [28]. A PDE is said to possess the Painleve 
property if its solutions are single- valued about the movable singularity mani- 
fold. To be more specific, if the singularity manifold is given by F(z±, Z2, ■ • • , z n ) = 
0, then a solution of the PDE must have an expansion of the form of eqn.(5). 
Substitution of this expansion into the PDE determines the positive value of 
r (from leading order analysis) and defines the recursion relation for the Uj. 

Weiss et al. truncated the expansion at the "constant term " level [28], i.e., 



u = 



F~ r Y.u,F\ 



(5) 



3=0 



U = UqF 



+ Ui 



F- r+1 + --- + u m _ 1 F- 1 + u, 



(6) 
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Substituting back into the PDE, one obtains an over-determined system of 
equations for F and Uj. The benefit of the singular manifold method is that 
this expansion for a nonlinear PDE contains a lot of information about the 
PDE. 

Most nonlinear non-integrable equations do not possess the Painleve property, 
i.e., they are not free of movable critical singularities [29,30,31,32]. For some 
equations it is still possible to obtain single-valued expansions by putting a 
constraint on the arbitrary function in the Painleve expansion. Such equa- 
tions are said to be 'partially integrable' systems, as presented by Hietarinta 
as a generalization of the Hirota bilinear formalism for non-integrable sys- 
tems [33]. He conjectured that all completely integrable PDEs can be put 
into a bilinear form. There are also non-integrable equations that can be put 
into the bilinear form and then the partial integrability is associated with 
the levels of integrability defined by the number of solitons that can be com- 
bined in an iV-soliton solution. Partial integrability means that the equation 
allows a restricted number of multisoliton solutions. Berloff and Howard sug- 
gested combining this treatment of the partial integrability with the use of the 
Painleve expansion, truncated before the constant term level, as a transform 
to reduce a non-integrable PDE to a multi-linear equation [26]. The Berloff- 
Howard method is a powerful tool for solving non-integrable equations. In 
this section, we give an example to show how to obtain solutions using this 
method. 

Firstly, we consider the generalized RCSH equation 

u t + au xx + (3u xxxx + 7« - 5u 3 = , (7) 



and show how to find exact solutions of non-integrable dissipative partial 
differential equations by using the Painleve test and the Hirota multi-linear 
method. 

We take the transform truncated at the term before the constant term: 

u = u F~ r + Ul F' r+1 + ■■■ + u^F" 1 . (8) 



Analysis of the leading order terms gives r = 2. By substituting this expansion 
(8) into the RCSH equation (7) and equating the coefficients of the highest 
powers of F to zero, we obtain expressions for uq , u\ in terms of F, and these 
lead to the transform 



^ ^F xx /l20/? d 2 
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This transform leads to an equation which is quadrilinear in F, meaning that 
each term contains a product of four functions involving F and its derivatives. 



7 F 2 F 2 - 2FF t F 2 + 6aF* + 2F 2 F t F xt - 1 F 3 F XX + F 2 F t F xx 



■12aFF 2 F xx + 3aF 2 F 2 x - 90^F 2 F 2 X + - F 3 F 3 



x r xx ^ovLr -T xx o\)jjr x r xx -f vu/jr -T xx r r xxt 

2 Z7 Z7 l 1 OH Q Z73 Z7 1 OH ft Z7 Z7 17 Z7 I 1 n /Q Z72 Z72 

XXX 



+4aF z F x F xxx + 120PF*F XXX - 120(3FF X F XX F XXX + 10/3F Z F^ 
-Q.F F xxxx 30f3FF x F xxxx -\- 15/? F F XX F XXXX -\- Q(3F F X F XXXXX 



PF Fxxxxxx . (10) 

Substituting F — 1 + e 2kx+2uJt into eq.(10) and equating the coefficients of 
different powers of e to zero, we arrive at the one-soliton solution: 

where 



'2^/2 1/3 




'12) 



5- / Q~ 

and we require a = — ^ — as a constraint on the equation parameters. This 
solution is reminiscent of the sech 2 -solution, representing an ultrashort soli- 
ton at the minimum-dispersion wavelength, taking into account the effects of 
fourth order dispersion [20,34,35] of the equation 



Actually, the simple ansatz ip = u(x) exp(— ryt) reduces this equation to the 
stationary real cubic Swift-Hohenberg equation (7). 

Now we consider the generalized RQSH equation in the form [36] 

u t + au xx + (3u xxxx + 7« - 5u 3 - r]u 5 = . (14) 



Analysis of the leading-order terms gives us r = 1. By substituting eq.(8) into 
the RQSH equation (14) and equating the coefficients of the highest powers 
of F to zero, one obtains expressions for u in terms of F that lead to the 
transform 
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This transform again leads to a quadrilinear equation in F: 



F+ F 3 - F+F F 2 + 3F F 3 - 58F F F 2 - 10 8 F F F 2 

^ tx - 1 x ± ~ h J± xxxxx 1 ufj± xxxx- 1 x ± lu r J xxx A xx 1 - 

+20(3F XXX F 2 F + aF xxx F 3 + 308F 2 X F X F - Q08F XX F 3 
-3aF xx F x F 2 - 2^5 F 3 F + 2aF 3 F + 1 F X F 3 = . (16) 

On substituting F = \-\- e 2kx + 2l ^ t i n t the quadrilinear equation (16) and equat- 
ing the coefficients of different powers of e to zero, we obtain four algebraic 
equations, but these algebraic equations don't have solutions. 

We suppose a different transform: 

» = §. (17) 



This leads to an equation which is pentalinear in F and G: 

F 4 G — F X F 3 G — [3F XXXX F 3 G + &(3F XXX F X F 2 G — A(3F XXX F 3 G X 
+Q8F 2 X F 2 G - 3Q8F XX F 2 FG + 2A(3F XX F X F 2 G X - Q8F XX F 3 G XX 
-aF xx F 3 G + 2A(3F^G - 2A8F 3 FG X + 128F 2 F 2 G XX 
+2aF 2 F 2 G - A8F X F 3 G XXX - 2aF x F 3 G x + F 4 G t + (3F 4 G XXXX 
+aF A G xx - SF 2 G 3 - ^G 5 = . (18) 

Substituting F — 1 + e 2kx + 2u)t and G = 2u e hx+UJt into equation (18) and 
equating the coefficients of different powers of e to zero, we get the one-soliton 
solution: 

u = Uq sech(/cz) , (19) 



where 




V V J 

Substituting F = 1 + e 2kx+2uJt and G = u (-l + e 2kx+2uJt ) into equation (18) 
and equating the coefficients of different powers of e to zero, we obtain the 
kink solution: 

u = uq t anh(A;x) , (22) 
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where 



k = ±uq 



( 



24/9 



7/ 



) 



4 




+ 47?7 



2r) 



, a = 



.(23) 



Generally speaking, the application of this technique to complex variable equa- 
tions is difficult because of the complex leading order. However, we will show 
that this method works in some cases such as the complex Swift-Hohenberg 
equation. 

2.2 Direct ansatz method 

It is well-known that solutions of many nonlinear differential equations can be 
expressed in terms of hyperbolic functions like tanh or sech. This fact motivates 
the direct method of solution in which the starting point is a suitable ansatz so 
that the p.d.e. is expressible as a polynomial in terms of tanh or sech functions. 
Clearly this method is not as general as the Hirota method, since it will not 
work if there is no solution of the assumed form. However, in principle, the 
method is more straightforward than the Hirota method, and it is more useful 
in obtaining periodic solutions (in terms of elliptic functions). 

In this section, we give an outline of the direct ansatz method and illustrate 
facets of its application by considering some straightforward examples. 

First, we consider the RCSH equation (7) and show how to get exact solutions 
of non-integrable dissipative partial differential equations by using the direct 
ansatz method. The analysis of the leading-order terms gives us r = 2, and 
so we suppose the solution u = w sech 2 (foe). We obtain the soliton solution 
by substituting this ansatz into the RCSH equation. This yields an algebraic 
equation in terms of tan/Vfunctions by using formulae for hyperbolic functions 
(see Appendix), and we then equate the coefficients of different powers of 
the tanh function to zero. The solution is identical to that found by using the 
Hirota method. We find that the constraint needed on the equation parameters 
is 7 = T^ffl- Assuming this is satisfied, the solution is: 




(24) 



where 



a 



(25) 



20/9 
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Now we consider the RQSH equation (14). Analysis of the leading-order terms 
gives us r = 1, and so we suppose the solution u = Uosech(kx). Substitution 
of this ansatz into the RQSH equation yields an algebraic equation in terms 
of the tank function. Again, by equating the coefficients of different powers of 
the tank function to zero, we obtain the soliton solution, and note that it is 
identical with the Hirota method solution. 

We find that the constraint needed on the equation parameters is now 

9a 2 35 2 , 2y/6aS , x 

+ 7+ — ± ^ = 0, (26) 



100/3 ' 50r] Ih^Wl 

With this satisfied, we have: 

u = u sech(fcc) , (27) 

where 

In the same way, we obtain a kink solution by the ansatz u t&nh(kx) 

u = uq tanh(A;x) , (29) 



where 

k - ± ym> u °- 2 V ± 2 V ' m 

and the constraint needed on the equation parameters is 



3^7 + 2Vmn4 

OL= k . (oil 



It is easy to develop this method to get solutions in terms of elliptic func- 
tions. If we substitute a Jacobi elliptic function, as an ansatz, into the RQSH 
equation (14), we obtain an algebraic equation in terms of some Jacobi elliptic 
function. By using the formulae in the Appendix, and equating the coefficients 
of different powers of the elliptic function to zero, we then get the elliptic so- 
lutions. 

The elliptic function solutions of the RQSH equation (14) are as follows: 
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Jacobi cn function 



u — Uq cn(kx, q) , 



where 



k = ± 



a - 2aq ± ^(2? - l) 2 a 2 - Af3^A 



2(3A 



A = lQq 2 - 16q + 1 , 



i/, n 



2Pg(a + 10/3(2g - 1)A; 2 ) 



and the constraint needed on the equation parameters is 
24(3q 2 k 4 



V = 



UA 



Jacobi sn function 



u = Uq sn(A;x, q) , 



where 



k = ± 



a 



+ aq± l) 2 a 2 - 4/fy4 



2/L4 



A = g 2 + 14g + 1 , 



2k 2 q(a- 10/%+ 1)A; 2 ) 



and the constraint needed on the equation parameters is 
24(3q 2 k 4 



V = 



UA 



Jacobi dn function 



u = u dn(kx, q) , 
where 



k = ± 



(q - 2)a ± ^(q - 2) 2 a 2 - 4/3-fA 



A = q 2 - 16q + 16 , 
2 2k 2 (-a + 10(q-2)(3k 2 ) 



u 







(43) 
(44) 
(45) 



and the constraint needed on the equation parameters is 
24/5 A; 4 



V 



(46) 



In principle, these Jacobi elliptic function solutions can be also obtained from 
multi-linear forms, because each Jacobi elliptic function can be expressed by 
theta functions which are Hirota r-functions in multi-linear form [37,38]. How- 
ever, calculation by this method is tedious, so we do not provide details of it 
here. 

Next we consider an elliptic function solution of the RCSH equation. We know 
from Painleve analysis that the RCSH equation has a double pole. Thus we 
look for elliptic functions having double poles. We suppose 

u = u + p(kx) . (47) 



Substituting this ansatz into the RCSH equation (7) yields an algebraic equa- 
tion in terms of a Weierstrass p function. By using formulae from the Ap- 
pendix, and equating the coefficients of different powers of the Weierstrass p 
function to zero, we obtain the following elliptic function solution: 

u = u Q + p{kx) , (48) 



where 



k = ±(-^Y , u = ±^=, (49) 
V120/V ' v/3055 V ; 

_ 2(10P 1 -a 2 ) 2VMa 2 -5Pj) 

92 ~ 3^ " ± ' ( } 

From its relation to the Weierstrass a function (see Appendix), we know that 
the Weierstrass a function is a r-function of the Hirota multi-linear form. 
Thus, we can construct this solution by the Hirota multi-linear method. How- 
ever, we do not give the detail of this approach here. 
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We can find another elliptic function solution of the RCSH equation, 

u = Uo + cn 2 (kx, q) , (51) 

where 



k = ±[-^- 2 , (52) 



UOfJq 2 
6 + 19u + 15ul ± VA 



(53) 



8 + 38u + 60u 2 + 30ul ' 
A — —15uq - SOul - 3ul + 12u + 4 , (54) 

and the constraints needed on the equation parameters are 



a = 20/5(1 - (2 + 3u )q) , (55) 
7 = 8/3fc 4 (8 - (23 + 30u )q + (23 + 60w + 45u 2 )q 2 ) . (56) 



We note that all above hyperbolic-function solutions can be derived in the 
limit of q — > 1 of above elliptic function solutions. 



3 Painleve analysis of the complex Swift-Hohenberg equation 



It is difficult to use a full expansion because the CQSH equation possesses 
both a complex leading order and non-integer resonances, and it generally has 
no consistency conditions. To see this, the CQSH equation is first rewritten 
as the following system of equations: 



tyt + (Pi + ifhtyxx + (7l + il2)ll>xxxx + (A*l + V^MV 

+(^ + ^)WV = #, (57) 

+ (a - ip 2 )r xx + (ti - H2)r xxxx + - V 

+(!/!- w/ 2 )MV = -^*- (58) 

To leading order, we set if) and -0* as 

tj} ~ ^ ^~ ?1 , ^* ~ , (59) 
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where 6 and 6 are the leading-order exponents. Upon equating the exponents 
of the dominant balance terms in CQSH 



(71 + il2)ipxxxx + {v\ + iv2)\i>\ 4 ip ~ , (60) 
we find 

6 + 6 = 2. (61) 

To find 6 and 6 we must also equate the coefficients in front of these terms. 
From eq.(57) we get 

^l 4 = -jTS^ + 1 )te + 2 )te + 3 )(f) 4 . < 62 > 

and from eq.(58) 

l^ol 4 = -^^6(6 + 1)(6 + 2)(6 + 3) . (63) 

Eqs.(62) and (63) are combined to give 

^^6(6 + i)(6 + 2)(6 + 3) = ^^6(6 + 1)(6 + 2)(6 + 3) .(64) 

Eqs.(61) and (64) are now the two equations we need to solve for 6 and 6- 
The result is 

6 = 1 + ia , 6 = 1 - iot , (65) 
where a satisfies 

(a 4 - 35a 2 + 24) ( 7l z/ 2 - 72Z7) + 10a (5 - a 2 ) (7^ + 72 i/ 2 ) = 0. (66) 

(To have a = 0, we need 71^2 = 72^1 as the condition for a chirp-less solution.) 
The leading-order exponents 6 = 1 + ia and 6 = 1 — 2a are not integers unless 
a = 0. This means that the CQSH equation already fails the Painleve test at 
the first step if a 7^ 0. 

When 7^2 — 72^1 = we can have a = (or a = y/E). In fact this is the 
condition of the next section, giving the chirp-less solution. This plain solution 
is not possible for the quintic CGL equation. 
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The expansions for ip and ip* therefore take the form 



^ = (i/joF" 1 -^ + i) X F- ia + ■■■) exp(iKx + iQt) , (67) 
iP* = (ilj*F- 1+ia + ^F ia + ■■■) exp(-iKx - ittt) . (68) 

The resonances for the CQSH equation can be calculated from the recursion 
relations for the -0/s and i/jj's, but we do not include them here because they 
are not important for our purpose. 

We can expect the following dependent variable transformation by the above 
analysis, 

G G* 
^ = JT^ exp(«fit) , V* = jy— exp(-ittt) , (69) 



(In general, the transform should be if) = F u- ia exp(iKx + iflt), however, in 
our case, this transformation also gives solutions in the next section.) 

In the same way, we can obtain the dependent variable transformation of the 
CCSH equation, 

1> = ^ ex P(^) > ^ = exp(-int) . (70) 



In the next section, we show the existence of exact solutions by using these 
transformations. 



4 Exact solutions of the complex Swift-Hohenberg equation 

4-1 The complex quintic Swift-Hohenberg equation 

We consider the CQSH equation in [23], 

iipt + (/ + 2id)ip xx + idif) xxxx + (6 2 - ibi)\ip\ 2 ip + (-C2 + ici)\if)\ 4 if) 

= (— a 2 + i(ai - d))ip , (71) 

(In the numerical work of Sakaguchi and Brand, a<i was 0.) 

Substituting the transformation (69) into the CQSH equation (71), we obtain 
a pentalinear equation 
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(id - Qai + a 2 - i)F 4 G + (6 2 - ib 1 )F 2 G 2 G* - (c 2 - zcO^G* 2 
—iF 3 F t G + iF 4 G t + (2/ + 4id)F 2 F 2 G + 24idF 4 G 
-{2f + 4id)F 3 F x G x - 2AidFF 3 G x - (2id + f)F 3 F xx G 
-36idFF 2 F xx G + 2UdF 2 F x F xx G x + 6idF 2 F 2 x G + (/ + 2id)F 4 G xx 
+ 12idF 2 F 2 G xx - QidF 3 F xx G xx + 8idF 2 F x F xxx G - AidF 3 F xxx G x 
AidF F X G xxx idF F XXXX G -\- idF G xxxx — . (72) 



Putting F, G and G* as polynomials in terms of exp(kx + cot) and substituting 
these functions into eq.(72) and equating the coefficients of different powers 
of e to zero, we obtain the following exact solutions. We can also obtain same 
solutions by using direct ansatz method. 

Bright Soliton 

ifj = gsech(kx)e int , (73) 



where 

fc = =h 



\^—> Q -- f + ^7T + a2 ' lgl ~ VJh • (74) 



The coefficients must satisfy the following relations: 
- 5a/5oi) 66; 

7 ,Cl = "7 



, ,2(4d-5y/dad M 2 d 

h = & 2 - , ci = Jr , c 2 = . (75) 



Dark Soliton 

i) = gt&nh(kx)e int , (76) 



where 



j fb 1 + 2db 2 2 2 2 A 2 



The coefficients must satisfy the following relations: 

n jr 2 

fll = rf(16A; 4 - 4A: 2 + 1) , c ± = , c 2 = . (78) 



Chirped Bright Soliton 
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V> = gsech(kx)e- ialosisech{kx)) e in \ 
where 



fc = ± 



\ 



fa + d(a 2 - 1) ± sjd{d - ai )A + (/a + d{a 2 - l) f 



dA 

A = a 4 - 6a 2 + 1 , 

n = k 2 (f - fa 2 + Ada + Ada{\ - a 2 )k 2 )) + a 2 , 
2 _ A; 2 (2da(23 - 7a 2 )k 2 + 6da - f(a 2 - 2)) 
151 ~ b 2 ' 

The coefficients must satisfy the following relations: 



k 2 (3fa - 2d(2 - a 2 ) - 2(2(10 - 19a 2 + a 4 )*; 2 ) 
Oi — 
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ci = — b 
c 2 = -b 



\9\ 

d{a 4 - 35a 2 + 24) 



2 (2da(23 - 7a 2 )P + 6da - /(a 2 - 2)) 2 ' 
10da(a 2 - 5) 
(2da(23 - 7a 2 )k 2 + 6da - /(a 2 - 2)) 2 ' 



Chirped Dark Soliton 

^ = gt&nh(kx)e- talogiscch(kx)) e int , 

where 



, 2_ _ 3/a v /(4rf-3/a) 2 -4rf(d-q 1 )A 

' V A 2r/. 1 2r/,l 
A = 16- 15a 2 , 

Q, = k 2 ((-2f - 6da) + 30dak 2 ) + a 2 , 



\9\ 



2 k 2 (2f + 6da - /a 2 + 10da(a 2 - 8)k 2 ) 



The coefficients must satisfy the following relations: 
k 2 (-4d + 3/a + 2da 2 + 10(2(4 - 5a 2 )A; 2 ) 



\9\ 2 

dk 4 (a 4 - 35a 2 + 24) 10rfaA; 4 (a 2 - 5) 



c i = i m , c 2 



<?l 4 ' Z 
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Elliptic function solutions of the CQSH equation are the following: 
Jacobi cn function solution 

ip = gcn(kx, q)e lQ,t , 



(95) 



where 



k = ± 



d(l - 2q) ± sjd(-12dq(q - 1) + cuA) 



A = lQq 2 - lQq+ 1 , 

n = fk\2q-l) + a 2 ,\gf = ^. 

The coefficients must satisfy the following relations: 

24dq 2 k 4 



Adqk 2 (l + 5(2q-l)k 2 ) 

Ol — 7—T7, , Cl 

\g\ 



\9\ 



c 2 = 0. 



(96) 
(97) 
(98) 



(99) 



Jacobi sn function solution 

ip = gsn(kx, q)e %nt , 



(100) 



where 



k = ± 



d(q + 1) ± Jd(-12dq + a x A) 



dA 



A = q 2 + Uq + 1 , 



n=-fk 2 (q + l) + a 2 , \g\ 2 = 



2fqk 2 



The coefficients must satisfy the following relations: 
4dqk 2 {-l + 5{q + l)k 2 ) 2Adq 2 k A 

\g\ 2 \gr 



, c 2 = 



(101) 
(102) 
(103) 



(104) 



Jacobi dn function solution 

ip = gdn(kx, q)e lUt , 



(105) 



where 
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k = ±^ , (106) 

A = q 2 -I6q + 16, (107) 

2fk 2 

V = ~fk 2 (q -2) + a2, \g\ 2 = . (108) 

The coefficients must satisfy the following relations: 

4dk 2 (-l + 5(q-2)k 2 ) 24dk 4 
bi = rr- 9 , ci = — — ,c 2 = 0. (109) 

\gr \g\ 



4-2 the generalized complex quintic Swift- Hohenberg equation 



Now we consider the generalized complex quintic Swift-Hohenberg(GCQSH) 
equation 

iipt + Pipxx + l^xxxx + A#| 2 ^ + H^lV = iSif; , (110) 



where all coefficients f3 — f3\ + ifo, 7 = 7i + 272, A* = A*i + *A*2 an d z/ = z/i + 
are complex and 5 is real. 



This equation can be easily normalized by rescaling t' = fi-yt , x' = y J ir ^x, so 
that f3\ and /ii can be 1 if f3\ and /ii are non zero. 

Substituting the transformation (69) into the GCQSH equation (110), we ob- 
tain a pentalinear equation 

uG 3 G* 2 + (24 + 50m - 35a 2 - 10ia 3 + a 4 )^F 4 G 
+2(-6 - Ilia + 6a 2 + ia 3 )j(2FF 3 G x + 3FF 2 F XX G) 
+fiF 2 G 2 G* - 6(-2 - 3*a + a 2 ) 1 (2F 2 F x F xx G x + F 2 F 2 G XX ) 
-(-2 - 3*a + a 2 )([3F 2 F 2 G + 3^F 2 F 2 X G + A 1 F 2 F X F XXX G) 

+ a)(2(3F 3 F x G x + 6 1 F 3 F XX G XX + A 1 F 3 F XXX G X + 
4 7 F 3 F X G^ + iF 3 F t G + (3F 3 F XX G + 1 F 3 F XXXX G) 

- n)F 4 G + iF 4 G t + (3F 4 G XX + -fF 4 G xxxx . (Ill) 



Putting F, G and G* as polynomials in terms of exp(A;z + a;t) and substituting 
these functions into pentalinear equation (111) and equating the coefficients 
of different powers of e to zero, we get the following solutions. We can also 
obtain same solutions by using direct ansatz method. 
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First, we consider chirp-less(a = 0) solutions. From eq.(66) we need jiu 2 = 

72^1- 



Bright Soliton 

ip = g sech(kx)e int , 



(112) 



where 



k = ± 



-&±^/% + 4fry 2 
272 



Q, = k 2 (l + 7l A; 2 ) , |#| 2 = 2k 2 (l + IO71A; 2 ) . 
The coefficients must satisfy the following relations: 

24 7l A; 4 



2A: 2 (/? 2 + 10 72 A: 2 ) 

to = no ' "i 



and /3i = Hi = 1 (normalized coefficients). 
Dark Soliton 

tp = gtanh(kx)e int , 



2A l2 k 4 





, v 2 = 








\g\ 4 



(113) 
(114) 

(115) 



(116) 



where 



k = ± 



4\ 



± \fPl + 16^72 



72 



n = 2A; 2 (87iA; 2 - 1) , \g\ 2 = 2k 2 (2^ik 2 - 1) . 
The coefficients must satisfy the following relations: 

24 7l A; 4 



2k 2 {20j 2 k 2 - fo) 

^2 = rT- 2 , vi 

\9V 



\9\ 



2A l2 k A 



and Pi = iii = 1 (normalized coefficients). 
Chirped Bright Soliton 



(117) 
(118) 

(119) 



■ip = gsech(kx)e 



-ia log(sech(fca;)) ^iQt 



(120) 
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where 



k = ± 



\ 



2a + (a 2 - 1)(3 2 ± y/(2a + (a 2 - l)(3 2 ) 2 + ASA 

2A 



A = 4a(a 2 - 1)71 + (a 4 - 6a 2 + 1)72 , 
n = k 2 ((l-a 2 )+2a(3 2 ) 

+k\(l - 6a 2 + a 4 ) 7l + (4a - 4a 3 ) 72 ) , 
\g\ 2 = k 2 {{2-a 2 ) + 3ap 2 ) 

+2A; 4 ((10 - 19a 2 + a 4 ) 7l + (23a - 7a 3 ) 72 ) . 

The coefficients must satisfy the following relations: 



A*2 = 



k 2 ((2 - a 2 )(3 2 - 3a) 

M 2 

2/c 4 ((10 - 19a 2 + a 4 )7 2 - (23a + 7a 3 )7i) 



v 2 = 



A; 4 ((24 - 35a 2 + a 4 )7i - 10a(a 2 - 5)72) 
W 

A; 4 ((24 - 35a 2 + a 4 ) 72 + 10a(a 2 - 5)71) 



and Pi — iii — 1 (normalized coefficients). 
Chirped Dark Soliton 



i> = gtanh(kx)e- ialoe{scch{kx)) e tnt , 



where 



k = ± 



\ 



3a - 2/? 2 ± ^{3a - 2[3 2 ) 2 - 48A 



2A 



A = 30a7! + (15a 2 - 16)72 , 

tt = k 2 (-2 - 3af3 2 + (I671 - 15a 2 7i + 30a 72 )A; 2 ) , 

\g\ 2 

= k 2 {{a 2 - 2) - 3a[3 2 + 10(4 7l - 5a 2 7i + 8a 72 - a 3 7 2 )fc 2 ) 
The coefficients must satisfy the following relations: 



^2 
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v x = 



v 2 = 



k 2 ((a 2 - 2)p 2 + 3a + 10(47 2 - 5a 2 <y 2 - 8a^ + a 3 ^)k 2 ) 
A; 4 ((24 - 35a 2 + a 4 )7i - 10a(a 2 - 5)72) 
£; 4 ((24 - 35a 2 + a 4 )7 2 + 10a(a 2 - 5)71) 



and Pi = Hi = 1 (normalized coefficients) . 

Elliptic function solutions of the GCQSH equation are the following 
Jacobi cn function solution 

ijj = gcn(kx, q)e int , 

where 



_ A-2q)P 2 ±^l-2q) 2 P 2 + 45 l2 A 
2 l2 A 

n = k 2 {{2q-l)+ 1 iAk 2 ), 

A = 16q 2 - 16q + 1 , 

|^| 2 = 2gA; 2 (l + 107i(2 g -l)A; 2 ). 

The coefficients must satisfy the following relations: 



A*2 



2qk 2 {(3 2 + 10 72 (2g - l)k 2 ) 



\9\ 



24 7l g 2 A; 4 



Vi = 



2l A 



, v 2 



2A l2 q z k 



\9\ \9\ 
and Pi = /Ji = 1 (normalized coefficients). 

Jacobi sn function solution 

ijj = gsn(kx, q)e int , 



where 



i 



l + q)P 2 ±J(q + l) 2 P 2 + 46 l2 A 



fc = ±. 

2 l2 A 

Q = k 2 (-(q+l) +1 iAk 2 ), 
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A = q 2 + Uq + 1 , 

|#| 2 = 2gA; 2 (-l + 107 1 (g+l)A; 2 ). 

The coefficients must satisfy the following relations: 
2qk 2 (-f3 2 + l0 l2 (q+l)k 2 ) 

\g\ 

24 7l g 2 £; 4 24~f 2 q 2 k 4 

and Pi — fj,i — 1 (normalized coefficients). 
Jacobi dn function solution 

ijj = gdn(kx, q)e tQt , 

where 

(g - 2)02 ± y/(q - 2) 2 p 2 2 + 4S l2 A 
2 l2 A 

Q = k 2 (-(q-2)+ ll Ak 2 ), 

A = q 2 - 16q + 16 , 

|<7| 2 = 2£; 2 (l-10 7l (g-2)£; 2 ). 

The coefficients must satisfy the following relations: 

2k 2 ((3 2 -10 l2 (q-2)k 2 ) 
\9\ 

24 7l A; 4 24 72 A; 4 

W ,V2 ~ W ' 

/Ui = 1 (normalized coefficients). 
^.5 TTie complex cubic Swift- Hohenberg equation 

Now we consider the CCSH equation 

if> t = ([i + iv)ip + ifip xx - d(if> + 2ip xx + if; 
-(5 + ij)\ip\ 2 ip . 



k = ± 



\ 



A*2 = 
V\ = 

and Pi = 
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Substituting the transformation (70) into the CCSH equation (157), we obtain 
a hexalinear equation 



(-27 - 5)G 4 G* 2 - d(120 + I54ia - 71a 2 - Uia 3 + a 4 )F 4 G 2 
+2d(24 + 2Qia - 9a 2 - ia 3 ){AFF 3 GG x + 3FF 2 F XX G 2 ) 
+ {a 2 - hia - 6)(12dF 2 F 2 G 2 x + 12d{2F 2 F x F xx GG x + F 2 F 2 GG XX ) 
+{2dF 2 F 2 G 2 - ifF 2 F 2 G 2 + UF 2 F 2 XX G 2 + AdF 2 F X F XXX G 2 )) 
+(2 + ia)(l2dF 3 F xx G 2 x + 2AdF 3 F x G x G xx + 8dF 3 F x GG x 
—AifF F x GG x + \2dF F XX GG XX + UF A F XXX GG ' x ~\~ 8dF F X GG xxx 
+F 3 F t G 2 + 2dF 3 F xx G 2 - ifF 3 F xx G 2 + dF 3 F xxxx G 2 ) + 2ifF A G 2 x 
-{d-n-iv + iil)F A G 2 - AdF 4 G 2 x - QdF 4 G 2 xx - 8dF 4 G x G xxx 
-2F 4 GG t - 4dF 4 GG xx + 2ifF 4 GG xx - 2dF 4 GG xxxx . (158) 



Putting F, G and G* as polynomials in terms of exp(kx + ujt) and substituting 
these functions into hexalinear equation (158) and equating the coefficients 
of different powers of e to zero, we get the following solutions. We can also 
obtain same solutions by using direct ansatz method. 

Chirped Bright Soliton 

tj) = gsech. 2 {kx)e- alos(scchikx)) e im , (159) 

k = ± . 1 , \g\ 2 = -(120 - 71a 2 + a^dk 4 , (160) 
V« 2 — 10 

H = 4(9 + 8a 2 )dk 4 , / = —I2adk 2 , (161) 

„ 4/ 9N (154a - Ua 3 )dk 4 
n = v- adk 4 {96 - 12a 2 ) , 7 = ^ — - — . (162) 

\g\ 

Now we look for solutions of elliptic function. We suppose 

1P = (g + p(kx)y nt . (163) 



Substitution this ansatz into the CCSH equation (157) yields an algebraic 
equation in terms of a Weierstrass <p function by using formulas in Appendix, 
and equating the coefficients of different powers of a Weierstrass p function 
to zero, we obtain the following elliptic function solution: 

ij = (g+p(kx)y n \ (164) 



where 



24 



\120dJ ' y 
_ 4d(5/x - 3d) 

/ = 0, 7 = 0. 



±\ 



2d 



± 



155 ' 
4(d - 5//) 



f2 = f . 



-2d 
1553 ' 



(165) 

(166) 
(167) 



^.^ TTie generalized complex cubic Swift- Hohenberg equation 



We consider the generalized complex cubic Swift-Hohenberg(GCCSH) equa- 
tion 

iipt + Pi^xx + l^xxxx + A^IV'I 2 ^ = Sip , (168) 



where all coefficients (3 — (3\ + i(3 2 , 7 = 7i + « 7 2j A* = A*i + ^2 an d 5 = Si + «5 2 
are complex. 

This equation can be easily normalized by rescaling t' = [i x t , x' = y^px, so 
that f3i and /ii can be 1 if fli and fii are non zero. 

Substituting the transformation (70) into the GCCSH equation (168), we ob- 
tain a hexalinear equation 



fiG* 2 G 4 + (120 + 154ia - 71a 2 - Uta 3 + a 4 )^F 4 G 2 

+2(-24 - 26m + 9a 2 + ia 3 )-y(4FGF x 2 F x G x + 3FGF X 2 GF XX ) 
-(-6 - hia + a 2 )(12 1 F 2 F 2 x G 2 x + 2A 1 F 2 F X GG X F XX 
+12 1 F 2 F 2 X GG XX + (3F 2 F 2 G 2 + 3 1 F 2 F 2 XX G 2 + A 1 F 2 F X F XXX G 2 ) 
-i(-2i + a)F 3 (12 1 (G 2 x F xx + 2F X G X G XX ) 
+A{(3F X GG X + 3 1 F XX GG XX + 2 1 F XXX GG X + 2 1 F X GG XXX ) 
+iF t G 2 + [3F XX G 2 + !F XXXX G 2 ) + {-iS - Q)F 4 G 2 
+2{(3F 4 G 2 X + 3 7 F 4 GL + ^F 4 G X G XXX ) 

+2(iF 4 GG t + (3F 4 GG XX + -fF 4 GG xxxx ) . (169) 



Putting F, G and G* as polynomials in terms of exp(A;x + u;t) and substituting 
these functions into hexalinear equation (169) and equating the coefficients 
of different powers of e to zero, we get the following solutions. We can also 
obtain same solutions by using direct ansatz method. 

Bright Soliton 

i> = gsech 2 (kx)e int , (170) 
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where 



^±4^'"= — (171) 

M =6^ >M2 = A>7i = -_ >7iI = -_. (172) 

and Pi = = 1 (normalized coefficients). 
Chirped Bright Soliton 



i) = g sech 2 (kx)e~ ialosiscchikx)) e int , (173) 
where 



k = =h 



-2(a 2 + 2a + 10)(a 2 - 2a + 10)5 2 



6a(a 4 + 16a 2 + 96)A + (3 2 A 

A = a 6 + 10a 4 + 8a 2 - 640, (175) 
/c 2 (-(a 6 + 10a 4 + 8a 2 - 640) + 6a/3 2 (a 4 + 16a 2 + 96)) 

2(a 2 + 2a + 10)(a 2 -2a + 10) U { ] 

2 _ A; 2 ((-a 6 - 3a 4 + 94a 2 + 1200) + 4a(2a 4 + 33a 2 + 205)/3 2 ) 
l5 ' ~ 2(a 2 + 2a + 10)(a 2 -2a + 10) ' ( ' 

The coefficients must satisfy the following relations: 

k 2 ((-a 6 - 3a 4 + 94a 2 + 1200)/j 2 - 4a(2a 4 + 33a 2 + 205)) 

2|(7| 2 (a 2 + 2a + 10)(a 2 -2a + 10) ' { > 

71 = ( a2 - 1Q ) + 6 ^ 2 (179) 

n 2P(a 2 + 2a + 10)(a 2 -2a + 10) ' 1 ' 

<v = (a 2 - 10)& + 6a 

72 2P(a 2 + 2a + 10)(a 2 -2a + 10) ' 1 J 
and fli — iii — 1 (normalized coefficients). 

Now we look for solutions of elliptic function. We suppose 

iP = (g + p(kx))e int . (181) 



Substitution this ansatz into the GCCSH equation (168) yields an algebraic 
equation in terms of a Weierstrass p function by using formulas in Appendix, 
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and equating the coefficients of different powers of a Weierstrass p function 
to zero, we obtain the following elliptic function solution: 



*l> = (g + p(kx))e int , (182) 
where 



^120 7i y ,y v^30^' 20 7l ' 1 ; 

2f3 2 + 2072^2 ^ 2^30(1 + 5 7 i(^ + ^i)) , lfi .. 

92 = 5— = , 93 = ± t== • (184) 

372^2 45^-7? 



The coefficients must satisfy the following relations: 



^ = -I2O72A; 4 , P 2 2 = ~30^29 2 , (185) 
135 72 3 ^ 3 2 = -8PM + 5 72 5 2 ) 2 , (186) 

and fli — fj,i — 1 (normalized coefficients). 
We find another elliptic function solution, 

1 Jj = (g + cn 2 (kx,q))e in \ (187) 



where 



k = ± {m^-) • < 188 > 

Q — —Si - 8 7 iA; 4 (8 - (23 + 30g)q + (23 + QOg + 45g 2 )q 2 ) , (189) 

6 + 19g+15g 2 ±VA 

q 2(g + 1)(15^ + 150 + 4)' iyUj 

A = -15g 4 - 30g 3 - 3g 2 + 12g + 4 , (191) 

, = Lz«£ __!_ . ( 192 ) 

3g 607ig 
The coefficients must satisfy the following relations: 



H2 = -120 72 g 2 A; 4 , (193) 

& = 20 72 (1 -(2 + 3<7)<z), (194) 

5 2 = -8 72 A; 4 (8 - (23 + 30g)q + (23 + 60g + A5g 2 )q 2 ) , (195) 

and Pi — jii — l (normalized coefficients). 
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5 Conclusion 



A number of pattern formation phenomena are descrived by the (2+1) dimen- 
sional complex Swift-Hohenberg equation. Among them is an approximation 
of a Maxwell-Bloch system for the transverse dynamics of a two-level class 
B laser, filamentation in wide aperture semiconducter lasers etc. (1 + 1) di- 
mensional version of this equation is one of the basic equations for modelling 
temporal behaviour of laser systems. Among others, it describes, short pulse 
generaion by passively mode-locked lasers with complicated spectral filtering 
elements. The knowledge of its solutions presented in any form will help to 
understand better the processes in such systems. In this paper, using Painleve 
analysis, Hirota multi-linear method and direct ansatz technique, we studied 
analytic solutions of the (l+l)-dimensional complex cubic and quintic Swift- 
Hohenberg equations. We considered both, standard and generalized versions 
of these equations. We have found that a number of exact solutions exist to 
each of these equations provided that coefficients are bound by special rela- 
tions. The set of solutions include particular types of solitary wave solutions, 
hole (dark soliton) solutions and periodic solutions in terms of elliptic Jacobi 
functions and Weierstrass p function. Clearly these solutions represent only 
a small subset of large variety of possible solutions admitted by the com- 
plex cubic and quintic Swift-Hohenberg equations. Nevertheless, the solutions 
presented here are found for the first time and they might serve as seeding so- 
lutions for a wider class of localised structures which, no doubt, exist in these 
systems. We also hope that they will be useful in further numerical analysis 
of various solutions to the complex Swift-Hohenberg equation. 
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A Formulas of hyperbolic functions and elliptic functions 

Let Dn = dn(A;x, q), Sn = sn(kx, q), Cn = cn(A;x, q). 
Jacobi elliptic functions satisfy the following two relations: 
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Sn 2 + Cn 2 = l, (A.l) 

Dn 2 = 1 - gSn 2 . (A.2) 

We have the following differential formulas: 
and 

1 r>W"(<r\ 

4 + g-4(2 + 5g)Sn 2 + 24gSn 4 . (A.4) 



qk A Dn(z) 



Now the limit q — 1 reduces this to dn(fcr, f) = sech(kx) = S(x), sn(kx, 1) = 
tanh(A;x) = T(x). Then 

and 

-^ = 5-28T 2 + 24T 4 . (A.6) 



The other main function we need for the periodic solutions is Jacobi sn func- 
tion. We have the following differential formulas: 

fS^" 1 -'- 1 ' (A ' 7) 

and 

1 Q "" ( \ 

FW = 1 + 149 + s 2 -^ + DSn*M 

+24g 2 Sn 4 (x) . (A.8) 

Now the limit g = 1 reduces this to sn(kx, 1) = tanh(/cx) = Then 

FW- 2 <- 1 + r2 >- < A ' 9 > 
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and 

1 T""(x) 



k 4 T(x) 



1(2 - 5T 2 + 3T 4 ) . (A.10) 



We give the following formula of Jacobi cn function: 
1 Cn"(x) 



k 2 Cn(x) 

and 

1 Cn""(x) 



-l + 2gSn i (x), (A.ll) 



1 + Aq - Aq(2q + 5)Sn 2 (x) 



k 4 Cn(x) 

+2VSn 4 (x) . (A.12) 

Finally, we give the following differential formulas of Weierstrass p function: 

[p\z)f = 4[p(z)] 2 - g 2 p(z) - g 3 , (A.13) 
p"(z) = 6[p(z)} 2 - \g 2 , (A.14) 

p""{z) = 12 (l0[p(z)] 3 - \g 2 p(z) - g^ , (A.15) 

where g 2 and g 3 are constants. Weierstrass p function is connected to Weier- 
strass a function by 

p(z) = --j^log<7(a0. (A.16) 

Weierstrass p function can be expressed by Jacobi elliptic functions: 

/ \ ei - e 3 2 e 2 - e 3 

p(a;) = e 3 + -— , A; = d - e 3 , q = , (A.17) 

sn z (kx, q) t\ — e 3 



where 
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ei + e 2 + e 3 = , eie 2 + eie 3 + e 2 e 3 = -— , eie 2 e 3 = — . (A.18) 
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